Topological phase for entangled two-qubit states 
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We present an unambiguous characterization of the rotation group 5*0(3) biconnectedness topol- 
ogy using two-qubit maximally entangled states. We show how to generate cyclic evolutions of these 
states, which are in one-to-one correspondence to closed paths in 5*0(3). The difference between 
the well known two classes of such paths translates into the gain of a global phase of tt for one class 
and no phase change for the other. We propose a simple quantum optics interference experiment to 
demonstrate this topological phase shift. 
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Two-qubit states are the simplest quantum mechanical 
systems displaying entanglement. Maximally entangled 
states (MES) allow for a clear and measurable distinction 
between classical and quantum mechanical predictions, 
as exemplified by Bell type inequalities jlj. They also 
play a decisive role in quantum information processes, 
like teleportation or dense coding '2l|. In addition, it is 
known that any quantum computing protocol involving 
many qubits can be irnplemented by concatenation of 
one or two-qubit gates This is why much attention 
has been paid to two-qubit systems, their properties and 
characterizations. 

In the present work, we intend to show that two-qubit 
MES can also be used to display, even experimentally, 
the well known double-connectedness of the 50(3) rota- 
tion group. Since the latter is often evoked to explain 
the minus sign multiplying a spin 1/2 state after a 27r 
rotation, we shall first argue that this is to some respect 
an ambiguous statement. Indeed, take one qubit state 
\^[t^0)) = a|0) -I- /3|1) subject to the Hamiltonian 
H ^ (an effective magnetic field along the z axis). 

At time t, the state reads \^{t)) = e^^a |0) + P |1). 
Let us turn to the Bloch sphere representation of \'^{t)). 
The three coordinates are given by the expectation values 
of the Pauli matrices. As is well known, the expectation 
value of the spin precesses around the effective magnetic 
field. At t = 2tt/uj it has experienced a full 2tt rotation, 
and a tt phase for the wave function. The tt phase has 
been clearly demonstrated on several systems, starting 
with beautiful experiment on spin 1/2 neutrons On 
general grounds, the change of global phase 7 under a 
cycle decomposes into its dynamical part, 7^ (as derived 
from the time dependent Schrodinger equation) and its 
geometrical part, the Berry phase 7^, whose origin relies 
on the Hilbert space peculiar geometry. The latter phase, 
7g, equals half of the area^ bounded by the representing 
part of the Bloch sphere \^ . The dynamical phase jd is 
simply related to the time average of the Hamiltonian. 



In the above case of one qubit precessing in a magnetic 
field, we have 7^ = — ttcos^ and jg = — 7r(l — cos0), so 
the expected global phase ^ — Jd + Jg — — for an ini- 
tial state with |a| = cos 0/2 and |/3| = sin6'/2. Let us 
remark here that this tt phase is present even for 6 — 0, 
where it is of pure dynamical nature (with no preces- 
sion and therefore no rotation at all). We may therefore 
question the geometrical interpretation of the tt phase. 
It would be more correct to state that this phase has a 
geometrical and a dynamical component, being purely 
geometrical only when 9 = n/2, when the initial state 
is orthogonal to the magnetic field. However, even in 
that case, we find problems in relating the phase to the 
double connectedness of 50(3). Indeed, the latter prop- 
erty relates paths on the 50(3) manifold, stating that, 
under continuous deformations there are only two differ- 
ent classes of paths, forming a two-element Z2 first ho- 
motopy group: Hi(50(3)). Note that two dimensional 
rotations (50(2)), have a more complicated structure: 
Hi(50(2)) — Z JJ. The above experiment on a spin 
1/2 particle was done with a constant magnetic field and 
therefore a constant rotation axis. A natural question is 
the relation of the measured tt phase to topological prop- 
erties of the rotation group: is the tt phase related to the 
subtle nature of 50(3) or to a property shared by 50(3) 
and its subgroup 50(2), as it seems to be the case? In 
the latter case, this phase gain is therefore not directly 
related to the 50(3) double connectedness, but rather to 
a multi connectedness property shared by both groups. 

In the following we aim to develop a non ambiguous re- 
lation between the 50(3) topology and the global phase 
of a wave function. The idea is rather simple: we use the 
one-to-one map between 50(3) and the set of two qubit 
MES 1^ . A suitable set of Hamiltonians that keep con- 
stant the degree of entanglement is constructed, so that 
the two-qubit MES traces a path in 50(3) when the ori- 
entation of the Hamiltonians is changed. The evolution 
is such that the dynamical phase vanishes while its ge- 
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ometrical counterpart, in a sense that will be precised, 
also vanishes, except at points where it abruptly changes 
to TT. Note that such a tt phase, although not related to 
the SO{3) topology, has already been discussed in the 
context of entangled states 0- 

A pure general bipartite state can be written as 



1^-) = a |00) + p |01) + 7 |10) + (5 111) 



(1) 



where a, /3, 7 and S are complex coefficients satisfying 
the normalization condition. The complex concurrence C 
of this state is defined as C = 2{aS — Pj) and its norm 
equals the standard concurrence C for pure 2-qubits as 
defined by Woofers ^ . The concurrence C is a measure 
of the degree of entanglement for pure bipartite states, 
assuming the value C = 1 for maximally entangled states 
and C = for product states. For the case of MES, state 
can always be written (up to a global phase) in the 
general form 0: 



m = \l-{a |00) + (3 |01) - P* |10) + a* |11)) , (2) 

where the star denotes complex conjugation. The Hilbert 
space of all the MES can thus be defined as Hmes ~ 



{{a,(3)eCy\c 



1, and(a, (3) ^ (—a, —(3)} 



S /Z2 = 5*0(3), i.e., there is an one-to-one correspon- 
dance between SO {3) and the MES. This represents a 
sphere (in 4 d) whose opposite points are identified. By 
writing states in the form of Eq. ||2J), we can see that 
each MES can be denoted by a pair of complex numbers 
(a,/3). 

We will now study closed trajectories in 50(3), cor- 
responding to the time evolution of MES subjected to 
precise Hamiltonians. The net result is that, depending 
on the homotopy classes to which the path belongs, the 
MES will acquire either a or tt phase. In the following 
we will show how this purely topological characteristics 
can be theoretically implemented and illustrated via an 
interference experiment. 

Let's take as the initial MES the Bell state {a = 1,(3 = 

0) 



(3) 



I*> = \/2(|00> + |11»- 



In order to implement the cyclic evolution of state lO, 
we will consider simple polygonal trajectories in S'0(3): 
A sequence of evolution operators will act for a fixed 
finite time interval on the states, until they reach back 
the initial state Q. A geometrical picture illustrating 
this time evolution can be given by the hypercube of 16 
vertices (belonging to S^) depicted in Fig. To get a 
discretized approximation of 5*0(3) {— S^/Z2) we must 
identify opposite vertices, leaving altogether 8 distinct 
points, from A to H in Fig. and the edges connecting 
these points. They all correspond to MES. The bar over 



the letters denote the opposite points. Point A is our 
initial state ^ and point A is identified to —A. The 
first class of trajectories that will be studied starts at A 
and finishes at the same point, making a circuit in the 
polygon defined by points A^B^F^D^A in Fig. 
(Q. This class will be called from now on the plus class, 
since it is homotopic to the identity and does not imply 
in a global phase change of the initial state. An example 
of the other class of trajectory, the minus class, is given 
by the sequence of points A^B^F^E^A. It 
also starts at point A but ends at the point diametrically 
opposed to the initial one (A). Thus, there is a global 
phase of tt gained by the initial state. Physically, this 
phase depends on the parity of the number of times the 
state crosses the space orthogonal to the initial one. 

In order to keep the one-to-one relation between en- 
tangled states and 5*0(3) elements, we will consider, in 
our trajectories, evolution operators leaving the com- 
plex concurrence C constant. Taking huj = 1, such 
unitary evolutions can be realized by the Hamiltonian 
H ~ ^fi.a, whith a = cTxX + ayij + a^z, cTx, Oy and CTz 
and n — n^x + Uyij + rizZ is a normalized vector giv- 
ing the orientation of the effective magnetic field. This 
Hamiltonian represents a rotation acting on only one of 
the qubits. The evolution operator corresponding to this 
Hamiltonian has the form: 



U 



-iHt 



cos : 



lUz sm : 



cos -{- iriz sin * 



(4) 

where n± = n^, ± iny. Defining e = e^^^^/V2, we 
have, for the points appearing in the hypercube of Fig. 

-> (1,0), B ^ (e,e), C -> (e*,e*), D ^ (e,-e), 
E (£*,-£*), F ^ {i,0), G ^ (0,1), H ^ (0,i). By 
changing the interaction time t and the direction of ro- 
tation of our two- level system, i.e., vector n, we are able 
to reach, departing from point A, all points in 5*0(3), 
representing all possible MES. Notice that all rotations 
are performed on one particle only. In order to find the 
exact Hamiltonians that produce each class of trajectory, 
we will chose, for reasons that will become clear in the 
following, to perform the first two rotations in the first 
qubit and the last two in the second qubit. Thus, the 
total operator acting on our states can be written in the 
form U ®Id for the first half of the trajectory and Id®lJ 
for the second one. One can easily check that in both 
cases the dynamical phases vanish. As for the geometri- 
cal phase, it can also be defined for open paths, as pro- 
posed by Pancharatnam 3|. This phase also vanishes in 
the present case, except when the state crosses the space 
orthogonal to the initial state, where it abruptly changes 
by TT. 

The parameters in equation Q can be easily found 
by comparing this matrix to the necessary rotations to 
implement each part of the circuits. By doing so, we find 
that each part of the trajectory takes a time t = 27r/3 and 
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the only difference in the Hamiltonians is in the direction 
to which ft points. These orientations {nx,ny,nz) read: 





v/T73(-i,-i,-i) 


B F 


yT73(i,-i,-i) 


F ^ D 


V^("l,-l,l) 


D^A 




F ^E 


v^(l,-l,-l) 


E^A 


Vi73(i,i,-i) 



For both classes of trajectories, the Hamiltonians coin- 
cide from A to F, from where the vector n will take a 
different value for each class of trajectory. 

We will now address the question of measuring such an 
effect in an interference experiment. In order to test the 
topological phase, one needs a setup enabling the inter- 
ference of an entangled state after performing one of the 
trajectories above with a reference one. A simple way to 
do that consists of interfering entangled photon pairs, as 
depicted in Fig. (|2Jl. Polarization entangled photon pairs 
are currently produced with high fidelity and efhciency in 
non linear optical systems ^ . Such twin photons propa- 
gate in different space modes, which will be denoted as a 
and b. Using this notation, the MES state emerging from 
a non linear crystal is j^-) = ^/T/2{\HaVb) + \VaHb)). 
Photons in modes a and b are detected by detectors Da 
and Db, respectively. Before detection, photons in mode 
b enter a Mach-Zender interferometer (see Fig.Q). In 
the first 50% — 50% beam splitter (BS 1), mode b is com- 
bined to mode c, which is in the vacuum state , |0c). 
After that, the total state is: 

I*) =. i ilHaVbOc) + t \HaObVc) + \VaHbOc) + i \VaObHc)) , 

(5) 

which is exactly what is needed to interfere two entan- 
gled states: BS 1 produced two pairs of polarization 
MES, each one of them entangled to one arm of the 
interferometer. An experiment producing the refer- 
ence fringes consists in the introduction of a variable 
phase factor e*"^ in arm b of the interferometer. Thus, 
before entering the recombining BS (BS 2), the state 
of the system is, taking into account the dephasing, \'^) = 
i {e'^ \HaVbOc) + i \HaObVc) + e'* \VaHbQc) + i |KOfci?c)). 
The passage through BS 2 changes the state into 

I*) = \{{e'^ - mHaVbQc) + ^ \VaHbQc)) + 

+ l)(|i?aObK) + iKOb-ffc))). (6) 

Coincidence detections in detector Da and Db produce 
a counting rate P = 1/2|(1 — cos0)|. The trajectories 
one want to implement are realized by photons in modes 
a and c. The space of polarizations, described by the 



Poincare sphere, is equivalent to the Bloch sphere [lOj. 
For polarizations, all rotations in the Poincare sphere can 
be accomplished via wave plates properly orientated in- 
troducing a dephasing in one of the polarization axis. 
The most general rotation matrix can be produced by a 
sequence of three wave plates: the first one introducing a 
dephasing of i\) in the vertical polarization axis, followed 
by the second one which introduces a dephasing of (5 and 
makes an angle d with the vertical polarization. Finally, 
a third plate introducing a dephasing of — V' in the ver- 
tical polarization can be added. As a result, a rotation 
matrix is generated which has exactly the same form as 
Q with the following correspondence: the parameter 8 
plays the role of time while for the components of vec- 
tor n we have — — sin 20 cos "0, Uy — sin 20 sin -0 and 
Uz = cos 20. We can thus introduce a set of plates act- 
ing in mode a and performing the common part of the 
trajectories, i.e., from point A to F. The remaining of 
the plates are put in arm b of the interferometer and 
their orientation determines the topology of the trajec- 
tories. The total state, before reaching BS 2, is then 
I*) = i(e'^ |i/„140c) |i?aObl4> +e'^ iK^ffcO,) + 

«(—!)" |VaObi?c))- The value of n (odd or even), depends 
on the nature of the trajectory. After passing through 
BS 2, the whole state transforms into 

1^) = ^((e^^ - (-l)")(|i?aHOc) + \VaHb^,)) + 

z(e'^ + {-\T)(\Ha^bVc) + iKOfci/e))). (7) 

If we now make the same coincidence counts in detec- 
tors Da and Db^ we obtain the fringes P — 1/2| ((—!)" — 
cos (/>) I . If the trajectory is of the minus class, a dephasing 
of one half period will be seen in the interference pattern. 
However, if it is of the plus type, the fringes remain at 
the same place. The choice of splitting the action of the 
rotations between the first and second qubit was made to 
stress the interest of using entangled states in the exper- 
iment. Since physics shall not depend on the order the 
detections were made, if we suppose that a detection in 
Da is made, the rest of the state can be seen as a statis- 
tical mixture of polarizations. This means that dealing 
with entangled polarization pairs or with a completely 
unpolarized beam would produce the same results in our 
experiment. However, the one-to-one map from photon 
states to 50(3) elements is only valid for MES and it is 
only in that case that we can claim probing a topological 
property of the rotation group. 

The gain of the tt phase factor, as mentioned before, is 
related to the number of times the MES crosses the space 
of its orthogonal states. An interesting variation of the 
experiment could well illustrate a different situation: one 
could think of trajectories which are not necessarily cyclic 
but, having crossed this orthogonal states space, already 
correspond to a gain of the tt phase. This would mean 
that the final states interfering in BS 2 will no longer be 
the same MES. As a consequence, there will be a decrease 
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in the fringe visibility, but the desired dephasing would 
still appear. 

Another interesting point to address is the case of 
states whose entanglement is not maximal. It can be 
shown that in that case, the tt phase also occur. How- 
ever, the relation between non MES and the SO{3) group 
is more complex |^, as well as the geometrical interpre- 
tation of their time evolution. 

The authors thank M. Brune, J. M. Raimond and S. 
Haroche for patient listening and fruitful comments and 
P. H. Souto Ribeiro for discussion of the experimental 
aspects. Laboratoire Kastler Brossel, Universite Pierre 
et Marie Curie and ENS, is associated with CNRS (UMR 
8552). 



* Electronic address: "Perola.Milman@lkb.ens.fr| 
^ Electronic address: mosseri@gps.jussieu.fr 
[1] J. S. Bell, Physics 1, 195 (1964); J. F. Clauser, M. A. 
Horne, A. Shimony and R. A. Holt, Phys. Rev. Lett. 23, 
880 (1969). For experiments, see A. Aspect, P. Grangier 
and G. Roger, Phys. Rev. Lett. 47, 460 (1981); G. Weihs, 
et al, Phys. Rev. Lett. 81, 5039 (1998). 



[2] M. A. Nielsen and L L. Chuang, Quantum Computation 
and Quantum Information, Cambridge University Press, 
Cambridge (2000). 

[3] H. Rauch et al, Phys. Lett. 54 A, 425 (1975); S. A. 
Werner et al, Phys. Rev. Lett. 35, 1053 (1975). 

[4] S. Pancharatnam, Proc. Indian Acad. Sci. A 44, 247 
(1956); M. V. Berry, Proc. R. Soc. London A 392, 45 
(1984); Y. Aharonov and J. S. Anandan, Phys. Rev. Lett. 
58, 1593 (1987). 

[5] C. Nash and S. Sen, Topology and Geometry for Physi- 
cists, Academic Press (1983). 

[6] R. Mosseri and R. Dandoloff, J. Phys. A: Math. Gen. 34, 
10243 (2001). 

[7] E. Sjoqvist, Phys. Rev. A, 62, 022109 (2000); M. Erics- 
son, et al, quant-ph/0206063 (2002). 

[8] W. K. Wootters, Phys. Rev. Lett. 80, 2245 (1998). 

[9] For some examples, see K. Mattle, et al, Phys. Rev. 
Lett. 76 4656 (1996); D. Bouwmeester, et al, Nature 
390 575 (1997); D. Bouwmeester, al, Phys. Rev. Lett. 
82 1345 (1999); M. Franga Santos, et al, Phys. Rev. A, 
64, 023804 (2002). 
[10] See, for example. Christian Brosseau, Fundamentals of 
Polarized Light, a Statistical Optics Approach, John Wi- 
ley ans Sons, Inc. (1998); J. J. Sakurai, Modern Quantum 
Mechanics, Addison- Wesley (1994). 




FIG. 1: Hypercube depicting the two topologically distinct 
trajectories. The plus circuit, corresponding to no phase 
change, is represented by the sequence of points A — > i? — > 
F ^ D ^ A. The minus circuit, corresponding to the gain 
of a TT phase, is the represented by the sequence A B ^ 
F ^E^A. 
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FIG. 2: Scheme for the experiment enabhng the measure- 
ment of the global topological phase. Polarization MES are 
produced in the non linear crystal in modes a and b. Photons 
in mode 6 are combined to mode c in a Mach-Zender interfer- 
ometer. A variable dephasing tj) can be introduced in mode b, 
while wave plates properly oriented (Pi and P2) act on arms 
a and c, respectively, performing the desired trajectories 



